Fractal Properties of Anomalous Diffusion in Intermittent Maps 
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An intermittent nonlinear map generating subdiffusion is investigated. Computer simulations 
show that the generalized diffusion coefficient of this map has a fractal, discontinuous dependence 
on control parameters. An amended continuous time random walk theory well approximates the 
coarse behavior of this quantity in terms of a continuous function. This theory also reproduces a 
full suppression of the strength of diffusion, which occurs at the dynamical transition from normal 
to anomalous diffusion. Similarly, the probability density function of this map exhibits a nontrivial 
fine structure while its coarse functional form is governed by a time fractional diffusion equation. 
A more detailed understanding of the irregular structure of the generalized diffusion coefficient is 
provided by an anomalous Taylor-Green-Kubo formula establishing a relation to de Rham-type 
fractal functions. 

PACS numbers: 05.45.Ac, 05.60.-k, 05.40.Fb 



I. INTRODUCTION 

Since more than a century normal diffusion, charac- 
terized by a linear increase in time of the mean square 
displacement (MSD) of an ensemble of moving particles, 
provides a paradigmatic example of a stochastic process. 
Here the MSD can be expressed by the ensemble average 
(x 2 ) ~ t@ with an exponent (3 = 1, where x holds for 
the position of a particle on the real line at time t > 0. 
However, exponents (3^1 are also possible yielding two 
important examples of anomalous diffusion: supcrdiffu- 
sion with j3 > 1 and subdiffusion with < (3 < 1. More 
recently, the importance of anomalous diffusive regimes 
was realized not only in physics but also in chemistry, 
biology and economics. These regimes were found in the- 
oretical models and experiments related, among others, 
to turbulence, amorphous semiconductors, porous me- 
dia, surface diffusion, glasses, granular matter, reaction- 
diffusion processes, plasmas and biological cell motility, 
see Refs. 0, i, 1 S, H, S, 0, H, H for reviews. Parallel 
to this development low dimensional deterministic maps 
attracted much attention as simple models of anomalous 
dynamics which can be understood anal ytic ally 110 . 1 111 . 
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A characteristic feature of the stochastic processes gen- 
erated by these maps is provided by the probability den- 
sity functions (PDFs) of the dynamical variables. In case 
of normal diffusion these PDFs exhibit Gaussian forms, 
whereas for subdiffusion they yield tails with stretched 
exponential decay and for superdiffusion Levy power laws 

m 

Several theoretical approaches have been worked out 
over the past few decades in order to explain anoma- 
lous diffusion. Perhaps the most famous one is continu- 
ous time random walk (CTRW) theory, pioneered by the 
work of Montroll, Weiss and Scher [23|. Their stochastic 
approach was later on adapted to sub- and superdiffu- 



sive deterministic maps [TT|, El, El, El, El, El H, H . 
Related maps were originally proposed by Pomeau and 
Manneville as simple models of intermittency [To| . A 
classification of their dynamics was provided particularly 
by Gaspard and Wang combining stochastic with dynam- 
ical systems theory [14(. Exploiting the deterministic 
properties of these maps, anomalous diffusion was further 
studied in the framework of the thermodynamic formal- 
ism [11 [2(| , by means of periodic orbit theory [ltil . 21.25] 
and by spectral decomposition techniques [23[ . Currently 
aging phenomena [H HI] , non-ergodic behavior (28| and 
infinite invariant measures [11] are in the focus of investi- 
gations demonstrating that these simple models provide 
ongoing inspiration for important new research. 

However, all of the above studies focused on specific 
values of control parameters only for which these maps 
are to a large extent analytically tractable. That gener- 
ally the dynamics is more intricate was shown by calcu- 
lating the parameter-dependent diffusion coefficient of a 
piecewise linear one-dimensional map, which turned out 
to be fractal [H, [H, HI] . Similar behavior was detected 
in more complicated models like the climbing sine map 
[HI], in experimentally accessible systems like particles 
bouncing on corrugated surfaces [341 ] and in models of 
Josephson junctions [351 ] . These findings can be explained 
by the existence of deterministic dynamical correlations 
that ar e top ologically unstable under parameter variation 
USE 111; for a review see Part 1 of Ref. 

In this paper we show that a fractal parameter depen- 
dence of physical quantities is not only typical for low- 
dimensional periodic deterministic dynamical systems 
exhibiting normal transport laws but also for anomalous 
dynamics. As an example we study parameter-dependent 
subdiffusion in a simple deterministic map, however, our 
main arguments hold for a broad class of anomalously 
diffusive systems. A specific feature of our analysis is 
that it employs a blending of techniques from stochastic 
theory and the theory of dynamical systems. That way 
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we wish to contribute to the microscopic foundations of 
a general theory of anomalous deterministic transport, 
which does not yet seem to be fully developed. 

Our paper is organized as follows: The model is intro- 
duced in Section [TTJ As a rate of diffusivity we choose 
a generalized diffusion coefficient (GDC), which gener- 
alizes the diffusion constant known from normal diffu- 
sion. In Section IIIII we present results from computer 
simulations showing that the GDC has a nonmonotonic, 
irregular dependence on control parameters of the map. 
We give a qualitative explanation of this phenomenon 
by arguing that the GDC is a self similar-like, fractal 
function. The GDC is furthermore conjectured to be ev- 
erywhere discontinuous. In Section llVI we briefly review 
CTRW theory. In Section [V] we show that an amended 
version of it correctly describes the coarse functional form 
of the GDC. Apart from fractal parameter dependencies, 
the map under consideration exhibits an interesting dy- 
namical transition, which is studied in detail in Section 
IVII In Section I VIII we compute the PDF of this map 
from simulations and analyze it by means of a time frac- 
tional diffusion equation. Section I Villi finally outlines 
a deterministic approach for analyzing the fractal GDC, 
which is based on a Taylor-Green-Kubo (TGK) formula 
for anomalous diffusion. This theory enables us to relate 
anomalous diffusion processes to fractal de Rham-type 
functions. Section IIXI summarizes our results. For an 
outline of this work we refer to Ref. 13811. 



II. THE MODEL 

We focus on a subdiffusive map which, restricted onto 
the unit interval, was introduced by Pomeau and Man- 
neville to describe intermittency [lfj |. 



x n+ i = M Zta (x n ) = x n + ax% 



< X n < - 



(1) 



where the parameter z > 1 holds for the degree of non- 
linearity and a > 1 is a second control parameter. In the 
following we will omit these two indices for convenience, 
M z . a = M. The variable i„el denotes the position of a 
point particle at discrete time n €E No- Translation sym- 
metry, M(x + 1) = M(x) + 1, and reflection symmetry, 
M(— x) — —M(x), complete the definition of the map 
on the real line. The iterated dynamics of this model 
is illustrated in Fig. [1] for particular values of the two 
control parameters: Given some initial condition xq, the 
equations of motion Eq. |T]) produce the next value x%, 
x\ determines X2, and so on. A typical trajectory of the 
system, displayed in Fig. [1] in a cobweb plot, is therefore 
represented in form of jumps xo — > x\ — > X2 — > . . . on the 
real line. 

For 1 < z < 2 the diffusion process generated by 
this map is normal, whereas for z > 2 it is anomalous 
[Til [Til This anomaly is due to the existence of 

marginally stable fixed points located at all integer val- 
ues of x n around which a particle gets 'trapped' for a long 




FIG. 1: Top: Illustration of the spatially continued map M 
Eq. |T} for parameters z — 3, a = 8 with a typical trajectory 
generated by the map (cobweb plot with dashed horizontal 
and vertical lines) starting at initial condition xq. Bottom: 
The demonstration of intermittency in a time series plot of 
a typical trajectory. Note the long laminar periods in the 
trajectory, which are interrupted by chaotic 'bursts'. Here 
and in the following figures all quantities are without units. 



time [15]. However, in regions with larger local deriva- 
tive, where the particle can 'jump' to different unit in- 
tervals, the dynamics becomes more irregular, see Fig. 
[1] Thus, a typical trajectory of the spatially continued 
map Eq. (TTJ consists of long laminar pieces interrupted 
by short chaotic bursts. Such a behavior is the hallmark 
of what is called intermittency jlQ |. An example of a 
typical intermittent trajectory generated by this map is 
shown at the bottom of Fig. [TJ 

In the groundbrea king works of Thomae, Geisel [ll| . 
Zumofen and Klafter [l7j | , main attention was paid to the 
time dependence of the MSD for specific fixed values of 
the parameter a. In contrast, our aim is to study how a 
suitably defined GDC behaves under variation of the two 
control parameters a and z. In Refs. [5, 39] an anomalous 
diffusion coefficient was introduced by 



D = D{z,a) = lim 



r(i + /3) (x 2 ) 



(2) 



FIG. 2: The generalized diffusion coefficient K, Eq. Q, for 
the map Eq. {JJ as a function of the map's nonlinearity z for 
different values of the prefactor a in the anomalous diffusive 
region z > 2. The semi- logarithmic scale is chosen in order 
to magnify the irregular fine structure. The values for a from 
bottom to top are a = 2.78, 3.14, 4, 5, 9.8. The curves consist 
of 100, 130, 200, 254, 450 points, respectively. 



where T is the gamma function and (3 is some constant 
that will be specified in the following. For convenience 
we define the GDC by 

2D 

K (^a) = — ; - . (3) 



K 



r(i + /3) 



For numerical simulations of the spatially extended map 
Eq. idJ) we have typically used an ensemble of 10 6 par- 
ticles, where the initial conditions are uniformly dis- 
tributed on the unit interval (0, 1). If not said otherwise, 
each trajectory was calculated up to 10 4 time steps. For 
different values of a we have confirmed that the numeri- 
cally computed power law dependence of the MSD is in 
full agreement with the CTRW solution [13] 



1, 



1 < z < 2 
z > 2 . 



(4) 



We remark that at the transition point between normal 
and anomalous diffusion a logarithmically corrected de- 
pendence of the MSD is obtained, 



(x\~n/ln.n for z 



(5) 



which again is in complete agreement with CTRW the- 
ory. Hence, for the rest of this paper j3 is denned by the 
theoretically correctly predicted values Eq. Q. 



III. NUMERICAL RESULTS 

We first fix the parameter a and study the dependence 
of K as a function of the degree of nonlinearity z. Nu- 
merical results are presented in Figs. and [3] if as a 



FIG. 3: The generalized diffusion coefficient if as a function 
of z for different values of a. The parameter values from 
bottom to top (at z = 2) are a = 3.14,4,5. The curves 
consist of 296, 360, 435 points, respectively. The inset shows 
several representative values of K around z — 2 for a = 5. 
Lines are guides for the eyes only. Different symbols from 
top to bottom correspond to simulations with 10 6 trajectories 
and the different computation times n = 10 2 (circles), 10 3 
(squares), 10 4 (rhombuses), 10 s (triangles) indicating ultra- 
slow convergence of K towards zero at z = 2. Also the curve 
calculated with 10 7 trajectories, each of the length of 10 4 
time steps, is shown. In the inset it coincides with the curve 
calculated with 10 6 trajectories, each of the length of 10 4 time 
steps. 



function of z appears to be highly irregular exhibiting a 
lot of structure on fine scales. Particularly Fig. [5] shows 
how the fine structure as a function of z evolves under 
variation of a. 

In Fig. [3] the calculations were performed both for the 
normal and for the anomalous diffusive regime. Eq. ([5]) 
suggests that at z = 2 the GDC exhibits a local minimum 
in form of a total suppression of the strength of anoma- 
lous diffusion, that is, K(2,a) — for all values of a. 
Interestingly, this transition from normal to anomalous 
diffusion appears to be approached continuously by K as 
a function of z. The reason is that for z ^ 2 logarithmic 
terms are still present in the dynamics, but they con- 
tribute for transient times only. This peculiar behavior 
does not only diminish the GDC but also significantly 
slows down the convergence of the simulation results. 
The inset of Fig. [3] shows several representative values 
of K in the transition region, where different symbols 
correspond to different computation times and different 
numbers of trajectories for fixed a. Due to the slow con- 
vergence, even for the largest computation times the re- 
sults for K are quantitatively still apart from the CTRW 
prediction, which holds in the limit of time to infinity. 
However, qualitatively there is a tendency towards zero 
at z = 2. This specific problem will be analyzed in full 
detail in Section [V] 

From normal diffusive maps it is known that iterations 



of the critical points of a map, which here are the points 
of discontinuity at a; = l/2+m,m 6 Z, play a crucial role 
in order to understand the complicated parameter depen- 
dence of a diffusion process [3(| HH . Accordingly, varia- 
tions of the height of the map h := M(i) = i + a(^) z 
strongly affect the parameter-dependent GDC, as we will 
discuss in detail below. Such variations are achieved both 
by changing z and a. In order to more clearly represent 
the impact of variations of z only on the GDC, we study 
if as a function of z for fixed height h. We emphasize 
that the topology of the orbit of the critical point is still 
affected by variations of z only, however, it is less sensi- 
tive to z than to varying h. Respective simulation results 
are shown in Fig. [U As expected, K(z, h) for fixed h is 
considerably smoother compared to Figs. [2] and [3J In- 
deed, the local minimum in the transition from normal to 
anomalous diffusion is now less obscured by a non-trivial 
fine structure. However, note that very small irregular- 
ities forming a self similar-like pattern are still present, 
as is suggested by Fig. @]and the inset, which depicts a 
magnification of a small characteristic region. The figure 
also shows the convergence of the data with respect to 
simulation time and number of trajectories from which 
the formation of the irregularities is clearly seen. 

Now we focus on the dependence of K on the param- 
eter h with the parameter of nonlinearity z kept fixed. 
In Fig. [5] the behavior of the GDC as a function of h 
is shown for different values of z which, as expected, is 
more non-monotonous than that of K(z). The structure 
of K(h) for different z looks somewhat similar, however, 
we could not detect a simple scaling. For z — 1 the map 
Eq. |T]) reduces to a piecewise linear map for which it 
was shown analytically and numerically that diffusion is 
normal and that the diffusion coefficient is a fractal func- 
tion of the slope a [3(| HH, [H, [H, H3, S3] , see the upper 
curve in Fig. [5J For z = 3 and over a larger range of a, 
K(a) as obtained from simulations is presented in Fig. 
[6] (a). Magnifications of the initial region, see Fig. [6] (b) 
and (c), show a self similar- like behavior of K on finer 
and finer scales indicating a fractal structure. 

A qualitative explanation for the fractality of all these 
curves can be given in terms of turnstile dynamics, which 
has successfully been applied in order to understand the 
fractality of the normal diffusion coefficient in piecewise 
linear maps. Here we outline the basic idea of this ap- 
proach only, for technical details we refer to Refs. [30l. l3ll] . 
The key ingredient of this method is the orbit generated 
by iterations of the critical point x c — 1/2 of the map M 
restricted onto the unit interval, x n = M n (x c ) mod I. 
If this orbit is periodic it defines a Markov partition 
with certain partition parts representing coupling regions 
(turnstiles) where, in the associated periodically contin- 
ued map, particles can jump from one unit interval to an- 
other one. Specific parameter values define Markov par- 
titions with specific turnstile couplings which, in more 
physical terms, generate specific sequences of forward- 
and backward scattering of particles that start nearby 
the critical point. It turns out that these Markov parti- 




FIG. 4: The generalized diffusion coefficient K as a function 
of z for the fixed height of the map h = \/3 as defined in 
the text. Different curves correspond to simulations with f 6 
trajectories and different computation times, from bottom to 
top on the right hand side of the figure n = f 2 , 10 3 , 10 4 . Note 
the inverse order of the curves at the point z — 2 marked 
by the vertical dashed line. Curves consist of 854, 646, 650 
points, respectively. Several values of K for n = 10 are 
also shown (circles). The data indicate clear convergence in 
time and the formation of irregularities. Additionally, several 
representative values of K calculated for 10 7 trajectories, each 
of the length of 10 4 time steps, are shown (squares). They 
coincide with corresponding values of K calculated with 10 6 
trajectories, each of the length of 10 4 time steps. The inset 
shows a magnified region of the curve calculated for n = 10 2 . 
The dashed box marks a region that appears to be self-similar 
to the whole region shown in the inset indicating the existence 
of an underlying fractal pattern. The curve was calculated for 
2 • 10 7 trajectories and consists of 200 points. 



tions are topologically highly unstable under parameter 
variation, in the present case both for varying z and a. 

In Fig. [5] (a) three parameter values are identified rep- 
resenting Markov partitions that correspond to two local 
minima and a local maximum of K(a) thus highlighting 
a specific fine structure in the whole curve. Higher-order 
Markov partitions can now be constructed (generated by 
higher-order iterations of the critical point) which are 
of the same type as the initial three identified in Fig. 
[5] (a), however, yielding new sets of parameter values. 
Such higher-order parameter triples are shown in Fig. [5] 
(b), (c). They identify the same type of fine structure on 
finer and finer scales. A similar analysis can be performed 
for other local peaks of K(a) qualitatively explaining the 
fractal structure of the GDC. More quantitatively, one 
may wish to calculate some fractal dimension for these 
curves. However, this is a hard task even if analytical 
formulas for a diffusion coefficient are available [4t| and 
is left for further studies. 

We conclude this section with an unexpected detail of 
the GDC. For the piecewise linear map at z = 1 the 
diffusion coefficient as a function of slope of the map is 
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FIG. 5: The generalized diffusion coefficient K as a function 
of parameter h for different values of z. From top to bottom 
it is z = 1,2.5,3,3.5,4. For clarity the curves are shifted 
vertically by 1.5, 1, 0.6, 0.4, 0.2, respectively. 



FIG. 6: (a) The generalized diffusion coefficient K over a large 
range of the parameter a at z = 3. Note that K(3, a) — for 
1 < a < 4. The dashed line represents the CTRW approxima- 
tion Eq. p8|l. The dashed-dotted line is the modified CTRW 
Eqs. (I31|l . (I34|l . (b) and (c) show magnifications of the initial 
region of (a) demonstrating the existence of irregularities on 
finer scales. The dashed line in (b) shows the modified CTRW 
Eqs. (|3U, (J33J). The sets of symbols included in (a) to (c) cor- 
respond to specific turnstile parameter values (see Section UTTl 
for explanations) . The inset in (a) depicts the development of 
a discontinuity in the GDC around a = 20. The different lines 
correspond to different iteration times, from bottom to top at 
left n = 10 3 (circles), 10 4 (squares), 10 5 (rhombuses), 10 6 (tri- 
angles), thus yielding a higher numerical precision. The curve 
with 10 7 trajectories and 10 4 time steps, which coincides with 
the one calculated with 10 6 trajectories and 10 4 time steps, 
shows that we gain convergence with respect to the number 
of trajectories. The data in (a) consists of 1198 points. 



continuous [321 ] . For z — 3, however, numerical results 
suggest that the dependence of K on a is discontinu- 
ous. In the inset of Fig. [5] we have magnified a small 
region of K around a = 20 (analogous behavior is found 
around a — 12, 28, 36, . . .). Note that for a > 20 particles 
can jump for the first time from the first unit interval 
to the fourth one. Increasing the numerical precision by 
increasing the computation time, K(a) near a = 20 grad- 
ually approaches a step function, which suggests that a 
discontinuity develops. Similar discontinuities have al- 
ready been reported for Lyapunov ex pon ents as functions 
of a bias in piecewise linear maps [41(. More compli- 
cated ones are known to exist for diffusion coefficients 
of nonlinear maps exhibiting bifurcation scenarios (33] ] 
and have recently been highlighted in research on trans- 
port in polygonal billiard channels [42| . An explanation 
why for our model the GDC is discontinuous for z = 3 
while it is continuous for z = 1 may be sought in the 
completely different character of the PDF of the map 
for z = 3, which develops non-integrable singularities at 
all marginally stable fixed points. Details of this pecu- 
liar behavior will be discussed later, see Fig. [5] and the 



analysis after Eq. (|52|) . Note that at suspected points of 
discontinuity of the GDC like a = 20 the critical point 
is getting mapped right onto these singular regions of 
the PDF at the very first iteration. Simultaneously, the 
orbit of the critical point exhibits a transition from for- 
ward to backward scattering under variation of a, which 
leads to a local maximum of the GDC in terms of turn- 
stile dynamics. For z = 3 this change in the microscopic 
scattering process is drastically amplified by the singular 
behavior of the PDF probably leading to a discontinuity 
in the GDC. In contrast, for z — 1 the PDF is not sin- 
gular but a non-differentiable step function [3l[ , thus the 
GDC is locally maximal but still continuous. It would be 
valuable to have a mathematical proof of this conjectured 
discontinuity phenomenon. 

The discontinuity at a — 20 is furthermore a boundary 
point of the specific fine structure identified in Fig. [SJ 
which was previously discussed in terms of turnstile dy- 
namics. As we have just emphasized, this parameter 
value is determined by a specific periodic orbit of the 
critical point of the map. However, as was outlined be- 
fore, one can now construct infinitely many higher-order 



6 



iterates of the same type of periodic orbits yielding new 
sets of parameter values. These parameter values are re- 
lated to the same type of structure and, hence, identify 
the same type of discontinuity on finer and finer scales. 
Furthermore, such parameter values are typically densely 
distributed on the parameter axis [3(1 [31j ■ This leads us 
to conjecture that K(a) shown in Fig. [5] is discontinu- 
ous on a dense set of parameter values (probably being 
of measure zero). Note that the argument carries over 
to K(z) in Fig. [5J where a careful study of the largest 
point of irregularity at z ~ 3.4 reveals the same type of 
discontinuity. 

So far we have focused on the fine structure of the GDC 
only. The next section proceeds with an understanding 
of its coarse functional form. As far as we know there 
is only one analytical result in the literature trying to 
predict the whole parameter dependence of the GDC of 
this map However, as we will show the respective 

calculation needs to be modified in order to match to 
computer simulation results. This necessitates to briefly 
review the whole approach by explaining our corrections. 



IV. CONTINUOUS TIME RANDOM WALK 
THEORY FOR MAPS 

The CTRW theory of Montroll, Weiss and Scher [27| 
has become a standard tool to model diffusion in in- 
termittent maps like Eq. (TTJ [ll|, This approach 
assumes that diffusion can be decomposed into two 
stochastic processes characterized by waiting times and 
jumps. Thus one has two sequences of independent iden- 
tically distributed random variables, namely a sequence 
of positive random waiting times Ti, T%, T3, . . . with PDF 
w(i), J w(t)dt = 1, and a sequence of random jumps 
Ci,C2,C3»--' with a PDF \(x), f^° oo X(x)dx = 1. For 
example, if a particle starts at point x — at time 
<o = and makes a jump of length at time t n — 
Ti+T 2 + - + T n , its position is x = for < t < T x = t x 
and x = (l + C2 + ••• + Cn f° r t n < t < t n+ \. The prob- 
ability that at least one jump is performed within the 
time interval [0,t) is then f Q dt'w(t') while the probabil- 
ity for no jump during this time interval reads ^(i) = 
1 — J Q dt'w(t r ). The master equation for the PDF P(x, t) 
to find a particle at position x and time t is then 



P(x,t) 



dx'X(x - x) / dt' w(t - t') P(x', t') + 



+V(t)8{x) 



(6) 



It has the following probabilistic meaning: The PDF to 
find a particle at position x at time t is equal to the PDF 
to find it at point x' at some previous time t' multiplied 
with the transition probability to get from {x' , t') to (x, t) 
integrated over all possible values of x' and t' . The second 
term accounts for the probability to remain at the initial 
position x = 0. It is easy to check that this equation 
yields a normalized PDF P(x,t). The most convenient 



representation of this equation is obtained in terms of the 
Fourier-Laplace transform of the PDF, 



P(k,s) = 



dx e ikx 



dt 



l P(x,t) , (7) 



where the hat stands for the Fourier transform and the 
tilde for the Laplace transform. Respectively, this func- 
tion obeys the Fourier-Laplace transform of Eq. ©, 
which is called the Montroll- Weiss equation [27| . 



P(k,s) 



1 — w(s) 



1 



s l-\(k)w(s) 



(8) 



The Laplace transform of the MSD can be readily ob- 
tained by differentiating the Fourier-Laplace transform 
of the PDF, 



(x 2 (s)) = / dxx 2 P{x,s) = - 



d 2 P(k,s) 



dk 2 



(9) 



fc=0 



In order to calculate the MSD within this theory, it thus 
suffices to know X(x) and w(t) generating the stochastic 
process. For one-dimensional maps of the type of Eq. (TTJ, 
using the symmetry of the map the waiting time distri- 
bution can be calculated from the approximation 



X n +1 



dx t 
~dt 



ax f 



x t < 1 



(10) 



where we have introduced the continuous time t > 0. 
Solving this equation for initial condition xq yields 



x t 



1 



— — a(z — l)t 



(11) 



We now define that a particle makes a jump when 
it leaves the unit interval. Note that this definition 
is different from the one of Ref. [ll|, which used the 
interval [— 5,5]. That the unit interval is the right 
choice for calculating the parameter-dependent random 
walk diffusion coefficient was shown for z = 1 in Refs. 
[I! [H H E3, From Eq. one can then obtain 
the time t that a particle spends on the unit interval be- 
fore making a jump according to the condition Xt = 1. 
The waiting time thus becomes a function of the injection 
point Xq, 



t(x ) 



1 



a(z- 1) 



1 



- 1 



(12) 



Accordingly, the waiting time PDF w(t) is related to the 
yet unknown PDF of injection points, 



W{t) ~ Pin(x ) 



dxo 



dt 



(13) 



Making the assumption that the PDF of injection points 
is uniform, P,„ ~ 1, the waiting time PDF is calculated 
to 



w(t) =a[l + a(z- l)t]~ 



(14) 
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for long enough times t, where normalization is used to 
obtain the prefactor. 

If jumps between neighboring cells only are taken into 
account, the jump PDF and its Fourier transform may 
be assumed as [17j 



X(x) = 6(\x\ - 1), A(fc)=cos(fc) 



Combining Eq. §§§ with Eq. (fT5|) leads to the Laplace 
transform [111 ] 



w(s) 



s(l — w(s)) 

For the following calculations it is useful to define 

1 



7 := 



z- 1 



z > 1 



(16) 



(17) 



Note that for z > 2, 7 is identical with (3 defined in 
Eq. dH). For 7 7^ 2 the Laplace transform of Eq. JTi|) 
reads 

10(a) = 7 (6s)' ) 'e bs r(-7, 6s) = 1 - (b S y e bs r{l - 7, bs) , 

(18) 

with 6 := j/a and the incomplete Gamma function 
T(a,x) :=f™dte-H a - 1 , 

For 7=1 the Laplace transform of the waiting time 
distribution is 



-I \ 1 2s 2. ,2s, 

w(s) = 1 e ■> Ai( — 

a a 



(19) 



where E\{x) is the exponential integral, E\(x) := 
-ei(-x) = J™ dt^. 

Here we are only interested in the long time behavior 
of the MSD, which corresponds to taking the limit s — > 
in Eqs. ([15]). fT9"|) with 7 being constant. For < 7 < 1, 
taking the lowest order of the expansion we get r(l — 
7, bs) ~ r(l — 7) and 



w(s) ~ 1 - (6s)T(l -7) 
Eq. HH) then yields 

1 



r 2\ _ 



1-7 



' ' r(i-7)6T 

and its inverse Laplace transform is 

, n> a 7 sin(7T7) „, 

(x ) = -= — iV^i 7 , < 7 < 1 

\ / 27T7 1 +')' 



(20) 



(21) 



(22) 



Similarly, for the normal diffusive case 7 > 1 we get from 
Eq. (fT5]) as the small s-asymptotics of the waiting time 
distribution 



w(s) ~ 1 



7-1 



(23) 



Using Eq. (|16p and the inverse Laplace transform yields 



<* 2 > 



7-1 



*, 7>1 



(24) 



In the remaining case 7 = 1, the expansion of the expo- 
nential integral is ei{—x) — \n(x) + 7, where 7 is Euler's 
constant. For small s Eq. (fl9|) reads 



w(s) ~ 1 + 6s ln(6s) 



(15) and Eq. (fT6|) gives 



ln(l/6s 



(25) 



(26) 



The inverse Laplace transform of the first term can be ob- 
tained by applying Karamata's Tauberian theorem [45j |. 
which relates the Laplace transform of a function given 
in the form /(s) - A^±±p-L(l/s) for s -> and 7 > -1 
to the inverse F(t) ~ At 1 L(t) for t — > 00. The function 
L must meet the requirement of being slowly varying, 
L ( ux ') — ^ l for « > as 1 -* 00. In our case we have 

1 



7=1, r(2) = 1 and the function L(l/s) — ^ nrgl) 
fies the above condition. For the MSD we thus get 



< >- 



at 



ln(at) 



satis- 



(27) 



In summary, we obtain for the GDC 



A' 



,7-1 



< 7 < 1 
7 = 1 
7 > 1 



(28) 



with 7 = -j-Lj-. For 7 = 1 we get K — because of the 
logarithmic term in Eq. (|27|) . 

Now we compare this result with our numerical simu- 
lations. As an example, A as a function of a determined 
by Eq. (|28p is shown in Fig. [S] (a), where 2 = 3. For 
z > 2, the difference to the result of Ref. [ll| is by a 
factor of two. This yields approximately the right slope 
for small a, however, it does not reproduce the trivial 
value of A (4) = when particles cannot escape from the 
unit interval anymore. Hence, further modifications are 
necessary. 



V. MODIFIED CTRW THEORY 

The easiest way to modify standard CTRW theory is 
by varying either the waiting time or the jump PDF. 
However, the waiting time distribution is straightfor- 
wardly determined by the model as a power law for all 
intermittent values of z. This suggests to change the 
jump PDF Eq. (|15p only, which in a first attempt one 
may write as 

X(x) = ^S(\x\-l) + (l-p)S(x) . (29) 

Here the second term reflects the fact that the particle 
can stay on the unit interval with probability (1 — p). By 
assuming that the density of particles is uniform on the 
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FIG. 7: K as a function of z for a — 5. The bold black line de- 
picts again the computer simulation results from Fig. [3] The 
thin line that is on top at z = 2 represents the simple ran- 
dom walk approximation Eq. (|36[) for i = 2. The dashed line 
corresponds to the modified CTRW approximation Eqs. (|32|l . 
(|34|) . The thin line that is in the middle at z — 2 displays the 
first term of the TGK formula Eq. (|53jl . The inset reveals the 
existence of logarithmic corrections in the MSD by plotting 
n 13 1 {x 2 ) with respect to Inn for a — 3.14 and different z close 
to the transition point at z = 2. The dashed lines correspond 
to z = 1.5, 1.9, 1.95 from bottom to top, the dashed-dotted 
lines to z = 2.5, 2.1, 2.05 from top to bottom. The thick solid 
line represents z = 2, the thin solid line is proportional to In n 
as a guide for the eye. 

unit interval, p is determined by the size of the escape 
region, p(a) = 2A(a) with A = | — x c , where x c is the 
solution of the equation x c + ax% = 1. For normal diffu- 
sion it is well-known that K rw — p(a) provides a random 
walk approximation for the diffusion coefficient which is 
asymptotically correct in case of nearest neighbor jumps, 
p(a) <C 1 }3Cil3ll.l36l.l43l|. For farther than nearest neigh- 
bor jumps this approximation is straightforwardly gen- 
eralized by weighting the integer jump distance squared 
with the respective probability, based on the escape re- 
gion, to perform such a jump |44j |. Assuming that one 
knows about the power law time dependence of the MSD 
Eq. ([2]) for this map, one can apply the same approxima- 
tion to anomalous diffusion. In Fig.[7]the result is shown 
for K(z) at a = 5. As one can see, this simple argument 
indeed roughly reproduces the coarse functional form of 
K(z) and yields the exact value at z = 1. 

Accordingly, the modification of standard CTRW the- 
ory Eq. (|29p works for nearest neighbor jumps only. We 
thus further amend the jump pdf Eq. (|29[) by introducing 
a typical jump length I, 

\(x) = Zs(\x\-l) + (l-p)5(x) , (30) 

with Fourier transform X(k) = p cos(lk)+l—p. This jump 
length we define either by the actual mean displacement 



or, under the simplifying assumption of integer jumps, 
by the coarse-grained displacement 

l 2 = {\[M(x)}\} , (32) 

where the square brackets hold for the largest integer 
less than the argument. Here the averages, denoted by 
curly brackets, are performed starting from a uniform 
initial ensemble of walkers with respect to the conditional 
probability density that particles leave the unit interval. 
The latter ensemble averages are then averaged over time. 
Repeating the calculations from the previous section with 
Eq. dUD, we obtain the MSD 

(£^H tl ; 0< 7 <1 

<* 2 > = ^ 2 i^y < ^ = 1 ( 33 ) 

{ V" at ' 7 > 1 . 

i = 1,2, where the typical jump length Li is given by 
Eqs. (|3T]>. ((32)) . Hence, our final CTRW expression for 
the GDC reads 



, 7=1 (34) 

, 7>1 • 

Fig. M (a) shows the CTRW approximation Eq. j34]) with 
jump length li, Eq. (|31[) . which describes the coarse de- 
pendence of the GDC very well over a large range of 
parameters. The CTRW result with integer jump length 
l 2 , Eq. (02) with Eq. (|3"I]). is shown in Fig.© (b). It gives 
an asymptotically exact approximation of the GDC for 
small values of a, however, it also holds well for larger pa- 
rameters. A quantitative comparison moreover indicates 
(within numerical accuracy) that this combination yields 
results which exactly correspond to our numerical find- 
ings for K at all integer values of h — M(h) > 0, which 
for z = 3 are equivalent to a = 4, 12, 20, .... This general- 
izes results obtained previously for the normal diffusion 
coefficient at z = 1 [U, |43| . 

We now focus on K as a function of z. Fig.[7]shows the 
modified CTRW approximation Eqs. (|32|) . (|34|) in com- 
parison with simulation results. In the strongly anoma- 
lous regime of large z the approximation reproduces the 
GDC from simulations very well, however, for smaller 
values of z there are obvious deviations. Let us first ex- 
plain the problem for z — > 1, the region around z = 2 
will be discussed in the following section. 

Let us recall that according to Eqs. ([SI),®, for z < 2 
the map exhibits normal diffusion in the long-time limit. 
One can further split this regime into 3/2 < z < 2, where 
the MSD shows transient anomalous diffusion which be- 
comes normal for long times while for 1 < z < 3/2 it 
represents purely normal dynamics [l4j. However, in the 
latter case it is well-known that the waiting time distri- 
bution is exponential [3lj |. 



h = {\M(x)-x\} 



(31) 



w (t) = je-^ 



(35) 
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where 7 is the escape rate, and not a power law such as 
the CTRW approximation Eq. (fl"4")) . In fact, in the limit 
of z — > 1 Eq. ([34)) leads to K = aplj, whereas the correct 
random walk result for normal diffusion reads [31, 43] 



K = pl\ 



(36) 



Note that for z — > 1 and by assuming a uniform density 
on the unit interval, in case of I2 <C 2 this equation yields 
K = p, whereas for / 2 > 2 we get K — l\ thus recover- 
ing the two simple random walk results mentioned ear- 
lier. Eq. (|36[) is indeed straightforwardly obtained from 
CTRW theory if one repeats the calculations with the ex- 
ponential waiting time distribution Eq. (|35|) and 7 = 1, 
i.e., by assuming that on average a particle leaves the 
unit cell after one time-discrete jump. We thus conclude 
that the CTRW result Eq. ((34)) gives only a reliable ap- 
proximation for the GDC if 7 < 2, whereas for 7 > 2 the 
simple random walk result for normal diffusion Eq. 
should be used. This is in line with Fig. [7] 



VI. SUPPRESSION OF THE GDC AT A 
DYNAMICAL TRANSITION 

Let us now understand the behavior of K(z) around 
z = 2, where the map exhibits a dynamical transi- 
tion [la , |47| from normal to anomalous diffusion, see 
Eqs. ([2j) , ((4j) . Interestingly, this transition is clearly visi- 
ble in the strength of the GDC, which to our knowledge 
has not been reported before. We now explain this phe- 
nomenon in detail. 

As mentioned in Section IIII1 around the transition 
point there are significant deviations between CTRW the- 
ory and the simulation results. These differences are dis- 
played in Fig. [7] At z = 2 the CTRW approximation 
forms a non-differentiable little wedge with K(2) = at 
the minimum, whereas the simulations yield K(2) > 0. 
We recall that by increasing the computation time there 
is some slow convergence of the simulation data towards 
the CTRW solution, see the inset of Fig. [31 but that we 
were not able to achieve quantitative agreement with the 
CTRW prediction of K(2) = 0. 

This very slow convergence of the simulation results as 
well as the logarithmic term in the MSD at z = 2, see 
Eq. ([3"3")) , suggest that the deviations between simulations 
and CTRW theory are due to logarithmic corrections in 
the MSD for parameter values around z — 2. Curiously, 
such terms are not present in Eq. ([3"3"| for z 7^ 2 (7 ^ 1). 
By refining our CTRW analysis, we will now show that 
such logarithmic terms indeed exist. However, for z 7^ 2 
they hold for large but finite times only, whereas for z = 2 
they persist in the limit of infinite time. As is obvious 
from Eqs. ([3"3")l , ([3~4"]) , for z = 2 the surviving logarithmic 
term leads to a full suppression of the GDC, whereas close 
to the transition point finite-time logarithmic corrections 
yield a gradual suppression of the strength of diffusion. 

We start again from the exact expression for the 
Laplace transform of the waiting time PDF w(s), 



Eq. (jT8j) . Combining this equation with the jump PDF 
Eq. ([3"0]) and Eqs. ([5]), © we obtain the Laplace trans- 
form of the MSD 



(z~2) _ p lfb-' r s- 1 -" i e bs T-' L (l - 7, bs) 



Pi* 



(37) 



We immediately drop the second term, since its inverse 
Laplace transform is simply constant. Let us focus now 
on the limits s — > and 7 — * 1. In case of < 7 < 1 the 
r function in Eq. (|37[) can be expanded to 

r(i- 7 ,6s) = r(i- 7 )-^L_^p( 



1 - 



7 ~[ (i + 1 - 7)j! 



(38) 

From this equation it follows that as 7 approaches 1 from 
below one may consider not only the first term as in 
Eq. ([20]) but also the second one, since both terms are 
divergent for 7 — > 1~. All other terms are nonsingu- 
lar and can be safely neglected. We further expand the 
second term of Eq. ([38]) to 



(bs) 



1-7 



1-7 



3 (l- 7 )ln(fc S ) 

1-7 



= — - + \n(bs) + £ 1 ~ 1} {Hbs)Y . 
1 3=2 J - 

(39) 

Note that here we have obtained a series of logarithmic 
terms which, as we shall see, provides the logarithmic 
corrections we are looking for. By imposing the condition 
that Kl — 7) ln(6s)| < 1, we may keep only the first two 
terms of Eq. 



(bs) 



1-7 



1 



]n(bs) 



1 — 7 1 — 7 
Finally, we expand the first term of Eq. (|38[) to 



(40) 



r(i- 7 ) = 



1 



1-7 



-7+E 



j=2 



(1-7) J 



dt e^'(lnt) J 



( 41 ) 

where 7 is again Euler's constant. Because we are in- 
terested in the small 7 limit, we keep only the first 
term of this expansion. In summary, we obtain for the 
Laplace transform of the MSD Eq. ([37)) in the limits of 
|(1 - 7)ln(6s)| < 1 and 7 -> 1" 



1 



(42) 



Inversion of the Laplace transform yields as the final re- 
sult 



<- 2 > 



pl\V 



&iT(l +7)ln(i/Z»)' 



t < t r 



(43) 



with t cr := 6e 1/(1_7) . For t > t cr and < 7 < 1 one 
can drop the second term in the T function expansion 
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Eq. flSE), and the asymptotic CTRW result Eq. ([33]) is 
recovered. Interestingly, t cr diverges when 7 — > 1 _ , and 
at 7 = 1 one thus arrives at the asymptotic i/lnt depen- 
dence of Eq. (021) ■ 

Analogous corrections are obtained for 1 < 7 near the 
dynamical transition point, 7 — ► 1 + . By applying the 
same arguments as above we get 



<* 2 > 



Plfr 



&tT(1 +7) ln(i/6)' 



t < t c 



(44) 



where i cr := fee 1 ^ 7 1 \ In the long time limit we again 
recover the CTRW asymptotics Eq. ([33]) g|. 

These analytical findings are in agreement with re- 
sults from computer simulations. In the inset of Fig. 
■nP I (x 2X j is shown as a function of Inn close to the transi- 
tion point at z — 2 for different values of z and a = 3.14. 
The logarithmic corrections are less obvious if z is suffi- 
ciently different from 2, where / (a; 2 ) quickly converges 
to 1/K. However, the logarithmic corrections are getting 
significant when z approaches the dynamical transition 
value. We have thus identified the precise dynamical ori- 
gin of the suppression of the GDC at z = 2. 

Dynamical transitions are quite ubiquitous in intermit- 
tent maps and have been widely discussed in the litera- 
ture in terms of the time dependence of the MSD. How- 
ever, it appears that so far no attention has been paid 
to a possible critical behavior of the associated GDC. 
Our results lead us to conjecture that suppression and 
enhancement of the GDC are typical for dynamical tran- 
sitions in anomalous dynamics altogether [lg|. 



VII. TIME-FRACTIONAL EQUATION FOR 
SUBDIFFUSION 



The great success of the CTRW approach is related to 
the fact that it not only predicts the power 7 but also the 
form of the coarse grained PDF P(x,t) of displacements 
[l7T |. Correspondingly the anomalous diffusion process 
generated by our model is not described by an ordinary 
diffusion equation but by a fractional generalization of 
it. Starting from the Montroll- Weiss equation and mak- 
ing use of the expressions for the jump and waiting time 
PDFs Eqs. IHl), ([151), we rewrite Eq. © in the long-time 
and -space asymptotic form 



s 1 P 



.7-1 - 



2c6 7 ' 



k 2 P 



(45) 



with c = r(l — 7) and b — 7/a. For the initial condition 
P(x, 0) = S(x) of the PDF we have P(k, 0) = 1. It is now 
helpful to recall the definition of the Caputo fractional 
derivative of a function G 



1 



r(i- 7 ) 
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and its Laplace transform [49j, [50j , 



dt 



.dr<G 



- s 7 G(s) - s 7 - 1 G(0) 



(46) 



(47) 



Noticing that the left part of Eq. (|45p precisely coincides 
with the Laplace transform of the Caputo derivative of 
the PDF and turning back to real space, we arrive at the 
time-fractional diffusion equation 



&iP{x,t) _ d 2 P 



dfy 



dx 2 



(48) 



with D given by the modified CTRW theory Eqs. ([3]), 
(|34p . Time- fractional diffusion equations of such a form 
have already been extensively studied by mathematicians 
[5l| . Note that in application to our model another ver- 
sion of such an equation was proposed in Ref. [24[ , which 
uses a Riemann-Liouville fractional derivative. It can be 
easily shown that both forms of time-fractional diffusion 
equations are equivalent under rather weak assumptions 
[52j . Yet two other forms of subdiffusive fractional equa- 
tions (however, not with applications to maps) were pro- 
posed in Refs. [53| and @. Again, after some recasting 
these equations yield Eq. ([48]) [54] . 

The solution of Eq. (|4"5)) is expressed in terms of an 
M-function of Wright type (49j . 



P{x,t) 



1 



-mU 



7 



(49) 



with £ := \x\j \f~Dt 1 1 2 ^ where we use the representation 
of the M-function 

1 00 < \i— 1 
M(^ 7 ) = -V^— -r( 7 i)sm( 7 m) . (50) 

7T * — ' II — 1 

i=l y ' 

This solution gives exactly the same asymptotics that 
was obtained in Ref. [l7| for small and large values of £, 



P(x,t) ~ t 



-7/2 



£^=4exp(-6i£^) , f>l 



(51) 

where a\, a? , b\ are some constants that are given explic- 
itly in Ref. gf. 

Similar to the modified CTRW approximation of the 
parameter-dependent GDC we now find that the analyt- 
ical PDF Eq. (|4"9"]l describes only the coarse scale of the 
PDF obtained from simulations but not its fine structure. 
This mismatch can be understood related to the fact that 
our model consists of a periodic lattice defined by specific 
elementary cells. Any elementary cell develops a char- 
acteristic "microscopic" PDF, characterized by the map 
Eq. (TTJ) restricted onto the unit interval, which exhibits 
singularities at the marginally stable fixed points. The 
inset of Fig. [5] schematically depicts such a PDF of an ele- 
mentary cell after a large number of iterations. Precisely 
this functional form yields the building block for the PDF 
of the periodically continued map, see the main part of 
Fig. [5] If one eliminates the fine structure by averaging 
over whole unit intervals, one obtains a coarse-grained 
PDF that is in excellent agreement with the analytical 
solution of the fractional diffusion equation Eq. (|4"5)) , see 
Fig. 1 
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FIG. 8: Main figure: The oscillatory structure shows the PDF 
obtained from simulations of the map Eq. |T} for z — 3 and 
a — 8. The simulated PDF was computed from an ensemble 
of 10 7 trajectories after n = 10 3 iterations. The crosses (x) 
represent the simulated PDF coarse-grained over unit inter- 
vals. The continuous line in the middle yields the analytical 
solution of the fractional diffusion equation Eq. (|49|l with the 
same z while D is taken from simulations. The upper and 
the lower curves correspond to fits with an M-function form, 
calculated with the same z but for different D, and a Gaus- 
sian form, respectively. More explanations are provided in 
the text, see Sec lVIII below. The inset depicts schematically 
the simulated PDF for the map Eq. Q modulo 1. 



Such an interplay between the "microscopic" PDF of 
a single scatterer and the "macroscopic" PDF of the spa- 
tially extended system was already reported for the nor- 
mal diffusive model at z = 1 I3U H3, [44| and was re- 
cently also studied in billiards [551 ]. However, in case of 
anomalous diffusion the PDF exhibits a further remark- 
able property: Fig. [5]shows that the oscillatory structure 
of the whole PDF is bounded by two different functions. 
The upper curve is of M-function type and connects all 
local maxima of the microscopic PDF. These maxima are 
situated in regions of the map where the dynamics is reg- 
ular, due to the marginally stable fixed points. The lower 
curve, on the other hand, is Gaussian and is determined 
by all local minima of the microscopic PDF. These min- 
ima are generated in regions of the map being far away 
from the marginally stable fixed points, where the dy- 
namics is locally most strongly chaotic [56[ . Fig. [5] thus 
nicely exemplifies the microscopic origin of anomalous 
diffusion in terms of intermittency. 



VIII. TAYLOR-GREEN-KUBO APPROACH 
AND FRACTAL FUNCTIONS 

We now turn back to the parameter dependence of 
the GDC. In Section [V] we have shown that an amended 
CTRW theory correctly describes the coarse dependence 
of the model's GDC, whereas the fractal fine structure is 



not captured by this approach. This reflects the fact that 
CTRW theory is a purely stochastic approach involving 
a randomness assumption between jumps, see the ap- 
proximations involved as outlined in Section IIVI and fVl 
whereas the origin of the fractality of the GDC lies in the 
existence of long-range deterministic dynamical correla- 
tions, see the analysis in Section HTT1 

This motivates us to propose an alternative approach 
for analyzing the GDC, which is based on the Taylor- 
Green-Kubo (TGK) formula [H HI for diffusion in 
maps. This theory has successfully been applied to the 
fractal diffusion coefficient of the normal diffusive map at 
z = 1, where it could partly be worked out analytically 
[3l| . to a nonlinear map with a complicated bifurcation 
scenario [33j j and to normal diffusion in billiards [3(| . An 
advantage of this analysis is that it systematically incor- 
porates dynamical correlations. A disadvantage is that, 
in contrast to CTRW theory, for the map under consid- 
eration we could implement it only numerically. 

The basic idea is to generalize the TGK formula for 
maps to anomalous diffusion. Following the usual deriva- 
tion [36|, [53|, the first step is to transform the MSD in 
Eq. ^ into sums over increments, or velocities, Vk '■= 

Xk+l ~ x k , 



< {x r . 



>=< Y v k Y vi > ( n ~* °°) ■ ( 52 ) 



fc=0 



For normal diffusion the second step requires to define 
the ensemble average expressed by the angular brackets 
in terms of the invariant PDF obtained for the map re- 
stricted onto the unit interval. In case of the map Eq. (Q]) 
this still works for z < 2, where the dynamics of the spa- 
tially extended map is normal diffusive. However, for 
z > 2 it is well-known that the map of the elementary 
cell does not possess a non-trivial invariant PDF any- 
more. Mathematical analysis [2^ shows that in this case 
there exist two physically relevant invariant measures, 
one which is concentrated on the marginally stable fixed 
points and one that lies in-between on the interval (0, 1). 
The first measure has still nice, so-called SRB proper- 
ties and yields a PDF in form of a ^-distribution on the 
marginally stable fixed points. The second one, however, 
is not normalizable anymore and hence is called an in- 
finite invariant measure. In other words, if one starts 
a computer simulation from an ensemble of points uni- 
formly distributed on the unit interval, the underlying 
stochastic process is not stationary and the PDF of an 
elementary cell, computed by a histogram method, does 
not converge in time to a well-defined invariant PDF 

EM- 

Consequently, in contrast to normal diffusion Eq. (|52[) 
cannot further be simplified by using time-translational 
invariance, and for the GDC Eq. © we have to stop at 



K = lim — 7 

n — >oc JiP 



n-1 



< Y v 'k > + 2 < vkVk + i > 



k=0 



k=Q 1=1 



(53) 
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In numerical simulations we find that the first term alone 
is already proportional to rfi . If the system is ergodic, as 
for z < 2 in our model, this term boils down to n < Vq >, 
and by neglecting any higher-order terms in Eq. (|53p 
we recover the random walk result for normal diffusion 
[3lL [3^ . 43] . For z > 2 the situation is again more compli- 
cated. Here only generalized ergodic theorems may hold 
[28l [29j . which is intimately related to the existence of 
infinite invariant measures as outlined above. Whether 
in this case a CTRW result such as Eq. (|3~4"|) can be ex- 
tracted from the first term is a non-trivial open question. 

However, considering the first term only as an approx- 
imation of K for all z, the numerical results are depicted 
in Fig. [5] The comparison of this approximation with 
the previous simulation results shows that this term pro- 
vides a first little step beyond a modified CTRW ap- 
proach, since it reproduces the major irregularities of K 
as a function of a and even follows the suspected disconti- 
nuities discussed in Section [Ml However, our numerical 
precision is not sufficient to conclude whether it yields 
exact values for K at integer heights. 

Eq. (|53p thus provides a suitable starting point for a 
systematic evaluation of the fractal structure of the GDC. 
Since the series expansion in Eq. (|53[) is exact, work- 
ing out further terms one should recover more and more 
structure in the GDC on fine scales [H, l3q . Interest- 
ingly, the second term in Eq. (|53[) can be understood 
as a velocity autocorrelation function that, in contrast 
normal diffusion, depends on two time scales. The exis- 
tence of a second time scale points to the phenomenon 
of aging in dynamical systems [24|, [26| . In our case this 
can be understood as a consequence of infinite invariant 
measures. 

Here we do not further pursue these questions but focus 
instead onto a direct link between the GDC and fractal 
functions as provided by the TGK formula, see Refs. I3ll 
I331 . 13H ] for normal diffusion. As was shown in Ref. [36j |. 
Eq. (f53"|) also holds if the velocities Vk are replaced by 
the integer velocities jk ■= [xk+i] — [xk]- It is now useful 
to start from K expressed again by Eq. (T5"2"|) , trivially 
rewritten as 



n— 1 n — 1 

K= lim -g < J2i k J2i l > 



(54) 



fc=0 



1=0 



The two sums suggest to define the jump function 

n 

Jn(x) -=^2jk , 

which satisfies the recursion relation [3lT ] 

Jn{x) = jo{x) + J n -i(M(x)) . 



(55) 



(56) 



The jump function gives the integer value of the displace- 
ment of a particle after n time steps, which started at 
some initial position x = xq. It thus contains essential 
information about the microscopic scattering process of 




FIG. 9: Comparison of the first term of the anomalous TGK 
formula Eq. (|53[) (lower line) with the computer simulation 
results from Fig. [13(a). 



a particle by showing how sensitively the displacement 
depends on initial conditions. Eqs. (f54|) . (|55p imply 



1 2 
K = lim — r < J„ 

n^oo Til 3 



t (x) > 



(57) 



Fig. [TO] displays the product of jump functions Jn( x ), 
which governs the anomalous diffusion process, for a rep- 
resentative parameter value. Clearly, as time evolves the 
structure of this function is getting more complicated. 
The GDC, in turn, is determined by the cumulative func- 
tion 



T 2 n {x) 



(58) 



where we fix the integration constants by the condition 
that T^(0) = and by requiring that the whole function 
is continuous. Integration of Eq. (|56[) would yield a recur- 
sive functional equation for T ? (x) , which is of the same 
type as the one derived in Ref. [33| • The solutions of such 
equations are generalized de Rham-, respectively gener- 
alized Takagi functions [3l|,[33[. However, at present such 
generalized de Rham-equations cannot be solved analyt- 
ically, hence we do not work out the details but instead 
compute T^(x) directly from simulations, according to 
the definition Eqs. ([55]) . (T55|) . Results for the first four 
iterations of T% (x) are shown in Fig. [TT] With respect to 
the construction of T% (x) , it may not come too much as a 
surprise that in the limit of n —> 00 this function exhibits 
a fractal structure. Interestingly, the GDC is obtained 
by integrating over this structure. If one starts from a 
uniform ensemble of initial conditions the result reads 



K 



lim -3^(1) 

71— i-OO fih' 



(59) 



This equation relates the GDC to a fractal function repre- 
senting the microscopic scattering process of our model. 
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FIG. 10: The first four iterations of the jump function Jn{x) 
defined by Eq. ([55} for z — 3, a = 8, n = 0, 1, 2, 3 based on 
a uniform ensemble of TV = 10 3 initial conditions. Note that 
there emerges a complicated fine structure. 
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FIG. 11: The first four iterations of the generalized de Rham- 
type function T%(x) as defined by Eq. (|58[) , which is obtained 
by numerically integrating J^ix) depicted in Fig. 1101 The 
parameter values are the same as in Fig. 1101 Note that there 
emerges a fractal structure. 



The numerical result for the GDC Eq. (|59)) is in good 
agreement with the one employing the MSD. In the light 
of Fig. [TT1 and Eq. ([55)1 , it may not be too surprising any- 
more that under parameter variation a highly non-trivial 
GDC is obtained from this model. 



IX. CONCLUSIONS 

In this paper we have studied subdiffusion generated 
by a paradigmatic one-dimensional intermittent map. In 
contrast to previous works here we focused onto the pa- 
rameter dependence of the GDC, which is an anoma- 
lously diffusive generalization of the diffusion coefficient 
known for normal diffusion. Computer simulations sug- 
gested that this GDC is a fractal function of control 



parameters. This finding was corroborated by a qual- 
itative explanation of the fractal structure in terms of 
complicated sequences of forward- and backward scatter- 
ing (turnstile dynamics), which are topologically unsta- 
ble under parameter variation. Our analysis furthermore 
led us to conjecture that the GDC is a discontinuous 
function of control parameters. 

In trying to understand the coarse functional form of 
the GDC, we applied standard CTRW theory to our 
model. By suitably amending previous calculations, we 
arrived at analytical approximations that enabled us to 
reproduce the whole coarse functional form of the GDC 
yielding asymptotically exact results in the limit of large 
and small parameter values. However, there are clear 
deviations between this theory and simulations around 
a dynamical transition from normal to anomalous diffu- 
sion. By refining our amended theory, we were able to 
explain these deviations in terms of logarithmic correc- 
tions leading to ultra-slow convergence of our simulation 
results and eventually yielding a full suppression of the 
GDC right at the transition point. These findings were 
confirmed by simulations of the MSD. 

We then studied in detail the PDF of our model. We 
first derived a time-fractional subdiffusion equation from 
CTRW theory. The coarse-grained PDF of our model 
turned out to be in excellent agreement with the non- 
Gaussian solution of our fractional diffusion equation. 
On a fine scale, however, the simulation results yielded 
an oscillatory structure reflecting the microscopic details 
of the intermittent scattering process. This structure was 
generated by the invariant density of the single scatterers 
hence revealing an interesting interplay between micro- 
scopic scattering and macroscopic diffusion. 

A more detailed understanding of the GDC was finally 
provided in terms of an anomalous TGK formula, which 
again is a generalization of the TGK formula for normal 
diffusion. The structure of this formula points at inti- 
mate relations to infinite invariant measures, aging, and 
generalized de Rham-type fractal functions. All these 
terms define very active topics of research. We thus hope 
that, along these lines, in future work more detailed rela- 
tions between mathematical infinite ergodic theory and 
the physics of anomalous diffusive processes can be es- 
tablished. 

In the present paper we have only studied a one- 
dimensional subdiffusive map, however, we expect these 
findings to be typical for spatially extended, low- 
dimensional, anomalous deterministic dynamical systems 
altogether. Further studies may also focus on the relation 
between the anomalous TGK formula and CTRW theory, 
on a spectral analysis of the Frobenius-Perron operator 
of this model, and on analyzing a superdiffusive map by 
similar methods. 
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